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f^ i Abstract 

(N 

oXj, In this paper, wc investigate the back-reaction of C/(l) gauge fields into a specific class 

'~, ' of inflationary settings. To be more precise, we employ a Bianchi-I geometry (taken as an 

anisotropic perturbation of a flat FRW model) within two types of Born-Infeld theories. 

^\ • Firstly, we consider pure Born-Infeld electromagnetism. For either a constant or a 6(0) 

coupling, inflationary trajectories are modifled but anisotropics increase; In particular, for 

fj ' the former coupling we find that a quadratic inflaton potential, within a constant ratio for 

Q^, the scalar and gauge energy densities, does not induce sufficient inflation, while in the latter 

the back-reaction in the cosmology determines (from the tensor-scalar ratio) a narrow range 

bJO! where inflation can occur. A Dirac-Born-Infeld framework is afterwards analysed in both 

non-relativistic and relativistic regimes. In the former, for different cases of the coupling 

p^ , (richer with respect to mere BI setups) between scalar and gauge sectors, we find that 

^ ' inflationary trajectories are modified, with anisotropy increasing or decreasing. In particular, 

0\ I a tachyonic solution is studied, allowing for a non standard ratio between scalar and gauge 

O^ ' matter densities, enhancing sufficient inflation, but with the anisotropy increasing. For the 

relativistic limit, inflationary trajectories are also modifled and anisotropics increase faster 

than in the non-relativistic limit. Finally we discuss how magnetic seed fields could evolve 

in these settings. 
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1 Introduction 



Modern cosmology has entered a new phase of precision, unthinkable only two decades ago 
[1]. However a fundamental issue remains open; namely the origin of galactic and cosmological 
magnetic fields. Such fields have coherence length of the order 1 Mpc, with a magnitude of 
around 10~^G [21 1^. The current explanation for the origin of such fields is that a small seed 
field emerged during primordial inflation, which was amplified at late times by the galactic 
dynamo mechanism [H [5] ; Cf . [6l O [HI [9l [10] and references therein for additional elements on 
cosmic magnetogenesis. As is well known from cosmological perturbation theory, the accelerated 
expansion stretches small scale physics to cosmological size, therefore it would appear natural to 
include a primordial magnetic field (e.g., a fluctuation) in models of early universe cosmology. 
A crucial feature in this setting is that a FRW cosmology with canonical fields is conformally 
invariant, therefore a mechanism to break the conformal invariance is necessary for magnetic 
fields to be present [11]. Proposals of solutions to this problem are discussed in |3l[5]. 

Additionally, dealing with a vector (gauge) field in cosmology may require some degree of 
anisotropjo |14j . Therefore we must go beyond the simplest FRW models to understand the 
dynamics of the universe with a magnetic field. In this context, let us mention the various 
anomalies present in the CMB, such as the suppression and asymmetries of the power spectrum 
in different hemispheres, together with the presence of a preferred direction due to the alignment 
of lowest multipole moments |15tll6j. These anomalies can be interpreted as a sign that isotropic 
statistics are only a leading order approximation, and that these anomalous measurements are 
due to non-isotropic expansion. Hence, it is reasonable to investigate (inflationary) anisotropic 
space-times with a non-zero gauge field. 

Model building with both magnetic seed fields and anisotropies has been explored in [17t \T8\ 
[T9\ [20l [2T] . One potential pitfall present in those references was that the magnetic field energy 
density grows rapidly during an inflationary phase, and was thought to destroy the inflationary 
dynamics. As pointed out in [18j , this claim was not based on concrete analysis and in fact 
was shown to be false. Inflation does occur in the presence of a gauge field and leads to weakly 
growing anisotropy - but is modified by the back-reaction on the inflaton |22[ 123] . This was 
due to the gauge field acting as a source term in the Klein-Gordon equation, whilst remaining 
subdominant to the scalar energy density in the Hubble parameter. The overall net effect was to 
further suppress the amplitude of magnetic field production after inflation. The above analysis 
refers to (canonical) Einstein-scalar-Maxwell theory. However within the context of open string 
theory, the effective action is non-linear and we can ask how the above results are modified by 
such non-linearities. 

Our purpose in this paper is to extend the discussion on gauge fields back-reacting on in- 
flationary cosmologies (and how magnetic seed fields can emerge and evolve consequently) to 
a class of non-linear theories described by Born-Infeld (BI) type Lagrangians. Within string 
theory, we have a Lagrangian description of the low energy regime of a (Dirichlet) Dp-hiane, 
which in a cosmological context gives rise to a class of theories known as Dirac-Born-Infeld 



^See ref. [121 113) for Einstein- Yang-Mills configurations consisting of FRW settings with non-Abefian (e.g.. 
SU{2)) gauge fields. 



(DBI) inflatioio. The inclusion of gauge fields in such models have been considered in |24^ [25] , 
but here we incorporate instead the ideas of anisotropic inflation and attempt to generalise the 
results. Moreover these papers assumed warped backgrounds which were asymptotically AdS, 
however more general classes of backgrounds are not asymptotically AdS which leads to non- 
trivial coupling between the gauge field and the dilaton. By assuming an extra condition, which 
we believe to be physically motivated, fixing the electromagnetic energy density in terms of the 
scalar energy density - we can consider more general configurations. 

Our paper is presented as follows. In section 2 we consider a simple BI theory, which 
reduces to the Einstein-Maxwell theory in a specific limit. In subsection 2.1 we take the usual 
BI configuration and in subsection 2.2 we analyse the gauge and scalar sectors through a b{4)) 
coupling. In both cases, we investigate how the gauge field back-reaction affects the inflationary 
dynamics. In particular, we describe how inflation can be modified and whether the anisotropics 
will increase or not. In section 3 we employ a phenomenological model of DBI inflation using 
a generalised expression for the DS-brane action. We consider the back-reaction of the gauge 
field in several different regimes. In particular, the energy density can be constrained so that 
inflation is driven by the scalar sector. In subsection 3.1 we discuss the non-relativistic limit, 
together with a tachyon configuration followed by a ^^5*5 solution. In subsection 3.2 we take 
the relativistic limit. In section 4 we conclude with a discussion of our results and suggestions 
for future work. In Appendix A we determine how magnetic seed fields can be generated and 
how DBI ingredients can contribute to the magneto cosmogenesis discussion. 



2 Born-Infeld theory 



In this section we consider the non-linear extension of the Einstein-Maxwell theory {U{1) 
gauge field) to the Einstein-Born-Infeld theory. Such a theory automatically includes electro- 
magnetic duality, and a causal structure that implies that the electromagnetic field is everywhere 
finite. Other works detailing non-linear Lagrangian densities include |26 [ I27 [ \28 \ I29j. A similar 
framework plays an important role in the effective description of D-branes within string theory, 
and therefore serves as a toy model for a more detailed analysis of open string dynamics. 

The matter Lagrangian for such a model can be written as follows 



^^,^^^^^-V{^) + l 




^ = -i^d^<l>d''(p - Vi(l)) + - ( 1 - ^/l + -F^pF-P ] , (2.1) 



where b is the BI parameter. One can see that as 6 — )• the theory is reduced so that it includes 
the Maxwell form; However the non-linear nature of the Lagrangian means that we can think 
of 6 as a deformation parameter which has mass dimension M~^. Cosmologically therefore, 
we expect b to deform the conformal nature of the gauge field term. We wish to study the 
implications of such a gauge field correction on the inflationary dynamics. 



''in the simplest of such models, our universe can be localised on a Z)3-brane which inflates as it evolves through 
a warped internal geometry. If such a model is to be correct, then there must also be a non-zero gauge field on 
the world-volume - corresponding to excited (Fundamental) F-strings. 



Using the residual gauge invariance, we consider solutions where ^o = and assume that 
the homogeneous vector potential A is aligned (and fixed) along the x-direction. This breaks 
the (spatial) S0{3) symmetry of the space-time to U{1) x SO{2), with the planar symmetry 
preserved transverse to the gauge field direction. Since our gauge field breaks isotropy, we must 
ensure that the background metric is anisotropic for the theory to be consistent [16^ \T8\ [30] . It 
is simplest to extend the FRW metric to a Bianchi I metric of the form 

where a{t) = e" is the scale factor of the homogeneous universe and a is herein a (perturbative) 
deviation from isotropy. Let us be more concrete about this point: In most of the case studies 
investigated in this paper, our starting point will be that of a spatially flat FRW space-time filled 
with a scalar field and a weak non-linear electromagnetic field; The energy density of the latter 
will be sought to be initially well below that of the scalar field matter to ensure that eventually 
the subsequent evolution is that of a Bianchi-I, but very close (in a perturbative sense) to a 
flat FRW. The above metric can thus be interpreted as describing a back-reacted geometry 
(sourced by the vector potential); When considering an initial state of isotropic inflation, one 
can consistently set A, o" — )• 0, in which case the metric reduces to that of flat FRW. 

The field equations for the gauge field derived from the Lagrangian can then be written 

= d^ (^V^Fatg^^'g^H^ + ^F^pF'^'r'^'^ , (2.3) 

however with the above ansatz for the gauge field, one sees that the Maxwell equation admits 
the following solution 

p2p-2a-8o- 
i2 _ f_Al (2 A) 

" " (1 + 6Ple-4"-4-) ' ^^■^' 

where Pa is an integration constant, which we can associate with the electromagnetic density. 
Note that in the limit 6 — )• the above solution reduces to that of Einstein-Maxwell theory. 
Because of the unique structure of the BI Lagrangian, there is a causal bound on the gauge field 
which can be written as 

bAj^ < e2°-4<^ (2.5) 

and if this bound is saturated, then the gauge field term completely drops out of the action 
leaving a residual cosmological constant term set by 1/6. In the isotropic scenario, the bound 
allows for a larger gauge field energy density at late times. Therefore the initial conditions for 
inflation are clearly sensitive to the magnitude of the gauge fleld. Calculation of the Einstein 
and scalar field equations gives us the following expressions, written in the on-shell formalism 
for the gauge field 

3«' - 3^' = ]j^ K<^' + ^ - r + r I ' (2-6) 




* + 3«^ = F77^|2-^^ ^ ' \, (2.7) 

2a + 3d2 + 3a2 = -i^ (^-±^2 ^ y + J- (^l _ ^l + foPJe-^-^.J J ^ (2. 

4> + 3a^ + V' = 0, (2.9) 



which form the basis of our analysis; Note however, that we assume b to be constant here. 
By relaxing this constraint the inflaton field equation becomes modified. This will be explored 
in section 2.2. In what follows, we will be interested in inflationary solutions of the Einstein 
equations in the regime of slow-roll and small anisotropics. 



2.1 Standard BI configuration 



Let us consider isotropic inflation as our initial setting. We therefore take the gauge field 
and a to zero, and the metric reduces to a flat FRW-form. Once the scalar potential is specified, 
the resulting dynamics can be solved for exactly. We will consider a small class of (chaotic) 
inflationary potentials for simplicity, such that V ~ mP(jf /p where p is unspecified. The scalar 
field equation can be solved in the Hamilton-Jacobi formalism to obtain 



02(a) = 00 - 2paM^ 



(2.10) 



Concretely for p = 2 we would then find that inflationary solutions using the WMAP normali- 
sation, correspond to m ~ 10~^Mp, where 4> ~ 0(10)^^,. 

Having established the isotropic inflationary trajectory, we can now ask how this is modified 
in the presence of a gauge field. Using the Einstein equations outlined above, we see that the 
equation for accelerated expansion (in its entirety) becomes 
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1- A/l + 6P^e-4«-4a 



(2.11) 



where inflation occurs when we neglect the kinetic terms compared to the potential. One can 
further see that inflation soon ends if the gauge field contribution dominates the scalar potential 
term in the slow roll regime - suggesting that the inclusion of a gauge field term will always 
spoil inflation. However in order to retrieve a more detailed appraisal, it is important to note 
that the gauge field is actually decoupled from the scalar sector, i.e., there is no source term 
in the equation of motion for - so the gauge field backreaction can only arise in the Klein- 
Gordon expression through the definition of the Hubble parameter. Let us then establish how 
a gauge field can modify inflation. For practical (computational) purposes, we will assume the 
slow roll approximation holds, and that initially a, o" ~ to leading order. Consequently, due to 
the special algebraic properties of the scalar sector, we can again employ the Hamilton-Jacobi 
formalism to solve for the field trajectory 
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b~\ 1 + bPJe-'^^ 



(2.12) 



valid for any scalar potential. Note that herein a prime denotes derivative with respect to a. Let 
us consider the case with p = 2, and consider a perturbative expansion |f| in Pa ■ The resulting 



""This is actually all we can do analytically, the full solution does however admit a numeric solution 



terms at leading order are therefore 

1 4aMp2 - (j)l 
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4' Ml 
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4' M^ 



+ ... 



" (2.13) 

where we have fixed the constant of integration to ensure that (^(0) ~ (^o and T are gamma 
functions |31] . The choice of sign arises from the square-root solution at leading order, although 
we will assume (for simplicity) that (f) is non-negative over the inflationary domain. Numerically 
we do not have to assume a perturbative expansion in the gauge field. However we see that 
the solution is remarkably insensitive to b. Indeed for fixed mass and field density, the scalar 
trajectories are almost identical even if b varies over factors 0(10^). 

Regarding the evolution of spatial anisotropy, under the herein assumptions, let us consider 
the second of the Einstein equations ()2.7p . in the limit where we can drop the a term. In the 
perturbative setting we have just discussed, we can re-write this equation for small o"^ as 






(2.14) 



a 3V{<t)) 

for an arbitrary potential. Note that this automatically implies that if V decreases during 
inflation, anisotropies are automatically increasing. Specialising to our solution for the quadratic 
case, we see that the back-reacted scalar field has higher order correction terms in P^ and 
therefore we see that the classical potential drives the anisotropies to increase. 

Because the gauge and scalar sectors are decoupled, the scale factor determines the sub- 
sequent evolution of electromagnetic energy. Indeed we can write the general energy density 
as 

PA-^l Ul + 6p2e-4a-4a _ l^ ^ (2.15) 

which suggests that at late times, the gauge field contribution to the energy density vanishes 
- although could still be important if it dominates the scalar contribution. In order to control 
this, one should arrange for the gauge field contribution to be sub-leading during inflation. We 
then consider the ratio of the two competing energy densities in the slow-rolling phase 



OA yl + bPle-^'^-^- -I 

Pct> 0^(0) 

which should be constant during inflation. Perturbative expansion of the ratio results in the 
general expression 

p2p— 4o-4cr 

which can only be constant when V{(1)) ~ Ce~'^°'~'^'^P\ due to the cancellation between powers 
of h. The constant C is dimensionless in the expression above. Clearly this solution does not 
correspond to the power law potential due to the relation between a and cj) through the Einstein 
equation. This means that we cannot simultaneously find inflating trajectories with a constant 
value of IZ for the power law solution. 



If we demand that the perturbative bound on TZ is the most important consideration, we must 
have V ~ Ce~^°P| in the isotropic hmit. The task now is to reconstruct the scalar potential 
as a function of <j). Note that this means the energy densities of scalar and gauge sectors are 
now proportional with C determining which is larger. Solutions with C > 1 imply that the 
scalar energy density is still dominant, whilst the gauge field dominates for C < 1. Referring to 
(j2.1ip , one sees that accelerated expansion of the universe requires C > 1 - therefore the scalar 
dominance should be considered the physical solution. Solving the Hubble expression for the 
scale factor yields; 

and then one can solve the scalar field equation to obtain 



r^ 



o + M,,/^ln(:^^] , A, = ^V1±^, (2.19) 



1 + C Vl + ^iV Mp 

which suggests the scalar potential must be of the form 



which is a decreasing function of (/> - because (p is initially chosen to be large. Note that this 
simplifies further if we assume the limit C >> 1. However this exponentially decaying potential 
is proportional to the gauge field term (which is already small), and moreover the functional 
form does not admit a slow-roll inflationary solution because the potential is far too steep. In 
fact, calculation of the slow-roll parameters give 

8(1 + C) 

which clearly do not allow for these inflating trajectories. This is not a surprising result. We have 
tried to constrain the theory to ensure that the total energy density remains constant during 
inflation - forcing the scalar and gauge field energy densities to be proportional to one another. 
As long as the scalar sector dominates, expansion occurs. However this constraint on the scalar 
potential reduces the Hubble friction in the equation of motion, thus the scalar field cannot 
drive a slow roll inflationary phase whilst satisfying the observational bounds. Thus although 
the universe expands, scalar field driven (slow roll) inflation is inconsistent. The simplest way to 
resolve this issue is to increase the number of scalar fields, ensuring that they all follow the same 
trajectory. This is known as Assisted Inflation [321 [331 [Ml E5] , whereby the combined effect of 
N scalars will be to increase the Hubble friction term in the equation of motion, ensuring that 
the overall centre of mass mode will follow an inflating trajectory. 

Imposing the constraint ()2.16p on the energy densities is too strong as it stands. A way to 
modify the theory is to introduce some kind of scalar source term, ensuring that this constraint 
is automatically satisfied. We will address this in the next section. 



2.2 6(0) Coupled solution 



Let us now consider a modification of tlie tlieory wliere we promote the BI parameter to be 
a function of the inflaton b{(p). This induces a non-trivial couphng to the gauge field, thereby 
allowing the gauge field to act as a source term as follows: 

<^ + 3d(A + d^V + ^ (l - ^1 + bPle-^'^-'^A . (2.22) 

For inflating trajectories the gauge field term eventually vanishes and we can solve the field 
equation once we specify the scalar potential as before. We will consider the quadratic potential 
as an explicit example. Because of the coupling to 5,^ ln(5) in the field equation, the backreaction 
on the scalar sector can occur before the gauge field energy density becomes comparable to 
the scalar energy density. Therefore we can solve the back-reaction problem at the level of 
the scalar field equation, and neglect the gauge field energy density in the Hubble parameter. 
Perturbatively we can then find the back-reacted Hamilton-Jacobi equation 

which is clearly sensitive to a term 5^ In b on the right hand side. Let us consider a few specific 
cases: 



• If (9(^ In 6 vanishes, or if b takes the form b ~ A(j)'^, then the source term simplifies dramat- 
ically. 

— In particular, there is a critical value of b for which the back-reaction term vanishes, 
and that is when b ~ 6o</>o/0^- This last condition is particularly interesting because 
it implies that at late times (when — t- 0) we recover the Einstein-Maxwell theory. 

— Inflationary trajectories with b ~ (j)~'^ are therefore unaffected by the additional 
electromagnetic energy density (for small Pa and with quadratic potential), however 
it may lead to an interesting contribution during reheating. 

• Consideration of the energy density ratio Tt = ^ (cf- ()2.16p ). i.e., meaning to find a 
constant ratio and where the energy density of the gauge field is sub-dominant, yields the 
following constraint upon the Born-Infeld parameter 

K'A) 1' Lple^f'-f'oVMl _ 2 (jl + b,Pl T V (2.24) 



(zfl + ^l + ftoPJ) 

where 6o is a constant set by the initial conditions on the inflaton field. Again we can only 
obtain a perturbative solution for the scalar field 
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(2.25) 



which is similar to the case of constant b discussed earher. This again suggests that the 
theory is insensitive to the Born-Infeld parameter. Numericahy we can integrate the full 
scalar field equation as before once we specify the scalar dependence of b{(f)). The resulting 
plots confirm that the solution is indeed insensitive to b{(p). 



Concerning the anisotropics generated during this regime, the expression is the same as that 
derived in (|2.14p and therefore we need not write it again. The scalar field solution contributes 
terms of higher order in Pa and thus we find the same result as in the decoupled case. 

The back-reaction on the scalar field is important cosmologically, however, because the num- 
ber of e-foldings is determined through the expression 



Let us assume that b ~ Acf)'^ for now, which is the simplest non-trivial case we can consider 
(aside from n = — 2). At leading order in a perturbative expansion we obtain the following 
solution for the scalar field at horizon crossing from ()2.23p 

cj)l ~ 2M2(1 + 2M) + ^e-<^°/^^' (l + ^) {Ei{2) - 7 - 2M^{1 + 2M)) + ..., (2.27) 

where (po is the initial value of the scalar field, and infiation ends ai <j) = \[2Mp and Ei is the 
expontential integral function [31]. We have also neglected higher order terms, and work in a 
regime where (j)^ /M^ » ln{(f)'^ / Mp) which allows us to neglect the logarithmic contributions. 
This latter approximation is physically well justified over the field domain. 

With this expansion we can estimate the cosmological perturbations at horizon crossing 
which we then write in terms of physical observables. It is convenient to define the following 
function (,{M) = 2 -|- 7 -|- 4J\f — Ei{2), which is positive definite for all M > 1. Then we find at 
leading order 

„ ^ 1 ! (2 + n)|(A0 /^e:^ 

^s - ^^2AA (1 + 2AA)2 (^ -^-2 I' i^-^«i 

r - '' [ 4(2 + n)g(AA) / ple-^g-^^r 

2 ^ il + 2Mfm' f (2 + <(AA)) ,-- 

' ~ 247r2M2 1 2(1 + 2AA) ' -2,^2 U' l^-^^i 



2m2(l + 2AA) L (2 + n)^(AA) 




n - 3^2M| y 4(1 + 2AA) ' -2, .2 ]}■ (2-31) 

Using the WMAP 7 year bound [35] on Ug, which satisfies rig ~ 0.963 it 0.012 |14|I36| and using 
AA = 60 to denote horizon crossing, we find that the following constraint on the backreaction 
parameter Q 

^ - m2M2 ~ -J2TW ^ ^ ^ 



where Q = corresponds to the case of zero gauge field. Now we note that the tensor-scalar 
ratio appears to be an increasing function of n, therefore the back-reaction term will lead to 
larger than observed values of r unless n is negative. This means that for successful inflation 
we can place a rough bound on the physically allowed values of n. We illustrate this in Figure 
1, shown below over a small range of n around n ~ —2 for the above bounds on the expansion 
parameter. As one can see, the lower bound implies a vanishing of the P^ term and therefore a 
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Figure 1: Plot of the scalar-tensor ratio as a function of power law index n. The WMAP bound 
implies that only values of n G (—2.4, —1.6) are physically acceptable. The lower bound arises 
because a physical r must be non-negative. 

constant value of r which is below the WMAP bound. This is to be expected, since this is just 
canonical inflation. For non-zero values of the Pa term up to the maximal, we see that there is 
a small window where r satisfies the WMAP bound, and is positive definite. 

Values of n outside the identified range are incompatible with observation, or non-physical 
therefore should be ruled out as viable candidates for infiation. As expected from our ear- 
lier observation of the Hamilton- Jacobi equation, the theory with n = — 2 is indistinguishable 
from infiation without gauge field (at leading order) and therefore easily satisfies the WMAP 
constraint. 



One could examine the entire parameter space to classify the infiating trajectories in terms 
of the BI functions, however this is not the main focus of this paper, and we leave the full 
parameter space evaluation for future work. We comment briefiy on the case of exponential 
coupling b{4>) ~ Ae"''^ , a is a constant, (cf. I2.24|) which may be of interest. In this case the 
scalar field admits the following perturbative solution 
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(2.33) 
where Ei{x) is the ExpIntegralEi(x) function |31j . and the field satisfies the usual boundary 



10 



condition that (/)(0) = (pQ. Note that this solution differs from that in (|2.25p through the use of 
the Ei function. We will now turn to a related issue, examining how gauge fields interact with 
the DBI theory of D3-branes. 



3 D-Brane theory 



As briefly mentioned in sections 1 and 2, the non-linearity type of structure present in the 
Born-Infeld frameworks turns out to provide the correct low energy description of Dp-branes in 
string theory. In this context a class of models within type IIB string cosmology, more often 
referred to as DBI-cosmology has been considered|j. Notwithstanding their interest, the simplest 
models of DBI-cosmology are plagued by several issues which disfavor them as inflationary 
candidates [371 [381 EH]. More complex models can be constructed which circumvent these issues, 
but they all appear sensitive to supergravity back-reaction in the relativistic limit [38l[l0l[ll]. 
All these particular models make use of internal spaces that are asymptotically AdS, and it 
has not yet been established whether more general backgrounds are more viable. Recently the 
non-trivial coupling between dilaton and gauge field has been exploited in [52] through the use of 
Wilson lines: Inflation in such a model appears to avoid several of the pathologies of previous DBI 
inflation, and suggests that such non-trivial couplings are important to realise an inflationary 
phase. Since such non-trivial couplings are generated through non-AdS backgrounds, our aim 
in what follows is to consider such geometries in the hope that they can be embedded into 
the full string theory. Throughout this section we will discuss how some simpler (Lagrangian) 
configurations can be of interest. We will show how we can compute modifications on the scalar 
field dynamics (and hence, the inflationary universe), given constraints on the electromagnetic 
energy density. More concretely, we will investigate specific limits of the equations of motion; 
Namely the non-relativistic and relativistic limit (cf. subsections 3.1 and 3.2, respectively). 
These will allow us to establish backreaction situations of gauge fields on the inflaton and 
possibly identify cosmological implications, namely observational. The feature to stress is that 
the allowed couplings induce herein a much richer class of dynamical possibilities than the 
content in section 2 (A brief discussion on generating magnetic seed fields suitable for these 
limits is first presented in Appendix A.). 

Let us start with the following effective theory for DBI cosmology for a single Z?3-brane 
(which is in Einstein frame), where we include a scalar potential which does not come from the 
world-volume calculation but could arise from a supergravity F-term, 

S = jd^x (-Fi{cP)^J-detiG^, + F3{<P)F^u) + qV^F^{<P) - V^Vi<P)] , (3.1) 

where the induced metric is 

G^, = g^u + F2{(l))df,(l)d,(l) (3.2) 



^Inflation in these models is driven by the motion of a probe D3-brane through a warped geometry. The 
warped geometry is typically taken to be a solution of the ten-dimensional field equations, which is then 'glued' 
to a compact manifold - and the probe brane is localised within the geometry far away from this gluing region. 
Since inflating trajectories are dependent on the particular warped background, we will consider instead a more 
phenomenological approach in this section - making reference to relevant string theory solutions as they arise 

11 



and g is a measure of the RR-charge carried by the brane, where q = 1 for a BPS D3-brane, 
q = —1 for a l)3-brane and q = corresponds to non-BPS branes such as employed in models 
of tachyon dynamics. We will assume that Fi is non- negative for simplicity. Negative tension 
objects must exist for a string compactification to be consistent, but we will not consider them in 
this particular paper. The bulk metric Qfj^u is assumed to be a solution of type IIB string theory 
in Einstein frame. Our ansatz ()2.2p is, in fact, general enough to extend to these backgrounds 
which typically have a non-trivial dilaton profile - which is important because the determinant 
contains additional powers of the dilaton, coupling to the Maxwell tensor in Einstein frame. 
Backgrounds falling into this class can be found in 



One can calculate the determinant at leading order and therefore can write the above action 
as follows 

S = J d^x (^-F,{cl>)V^ (^l + ^Fi{cP)F^,F^''^ + ,/^iqFiicP) -Vm^ , (3.3) 

where indices are raised and lowered with the full metric G^,y rather than g^,^. Note that we 
are neglecting any couplings to the axion, which in this theory could come in the form cqF A F. 
The Fi are functions of the scalar embedding, arising from placing the D3-brane in the non- 
trivial geometry. The factors of the brane tension are contained in the Fi parameter. This 
action appears superficially similar to that of Einstein-Maxwell theory, which typically has a 
coupling between scalar and gauge sectors as b^{4>)Ff^^F'^'^. In the D-brane theory we see that 
this coupling is fixed to be FiF^, which potentially allows for much richer solutions than the 
function b{(j)). Note however that the Maxwell tensors also couple to the inflaton through the 
term \/— G, which contain powers of (p^. For dynamic solutions, this coupling explicitly breaks 
conformal invariance. Therefore even for theories with FiF^ = 1, the conformal structure of the 
gauge field sector is still broken [25]. 



The modified Klein-Gordon equation is much more complicated, but takes the following form 



= Tl + ^F^.F^'A I FiF^^\l - 4e3") - 2F{ 



1 77 ^2 ^ e3-FiF2(^2 
1 - F2<p + 



:i - i^202) 



2\ll-F2(t>^ 



2FiF20e3" 1 + -^F.^F^'' - FiF^F^ i^^^ ^ ^ _ p^^2 ^ ^Sa^^^i 



+ qFJ-V'-'^ r^'^^'^^^JF^^Fn, (3.4) 



e3"FiF2F|(/> d 

4-^l-F2(/.2 

where a dot denotes a time derivative, and primes denote herein scalar field derivatives. For 
backgrounds which are asymptotically AdS, the string and Einstein frames coincide - leading to 
a theory with a trivial dilaton (at leading order). Canonical examples of such theories include; 

• Fi = F = (p^/X, Fa = /, F3 = 1/\/F, q = ±1 for the AdS solution; 

. Fi = y ((/)), F2 = 1, F3 = A, g = for the tachyon solution [IS Sg iZl HH] 
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We must now consider the Einstein equations for such a configuration. These will generally be 
complicated by the presence of the induced metric G rather than the space-time metric g. Using 
the variational identity 

SC^' = Sg^'^Sl + g'^Kj,)-^ - 5g^^g--K,{5>;, + g^^^K^r^Sl + 5'^" V)"'' (3-5) 

one can calculate the corresponding energy momentum tensor to be 



+ IfiFI^iTf^ [F^bFcdG'''{6'^ + g^'hec)-' + F^^F.^G'^^d'^ + g^^K,)-^ 



+ Fi^i - F2^^ [-G^bK + g'^Kb)-' + Gabg'^'K^iSj^ + g^^'h^brHs:^ + g'^hj^r'^ 

+ g^,u{qFi-V), (3.6) 

where we have defined hab = F2da4>db(j) for simplicity. Since the DBI action has a well-defined 
relativistic limit, we choose to define 7 = (1 — ^2^^)""^ as the analog of the 'gamma' factor in 
special relativity. The resulting Einstein equations arising from the Bianchi metric, can then be 
written as follows 

1 / FiF2 3 42p-2a+4CT\ 



-2°+4-Fi_F|A|7(72 + 1) 
6M2 



a + 3aa = ^tff '' (^.S) 



2a + 3a +3a = J^ [^ ' + V J (3.9) 

and one finds that the equation for an accelerating universe is given by 

a 2^ 1 / Fi(3-27g-7^) e-2"+^-FiF32^g7(37^ - 1) \ 

From the right hand side of this equation one sees that the energy density of the vector field 
is important when considering inflationary dynamics. Indeed we can identify a critical value of 
the gauge field which allows for acceleration (isotropic limit): 

Note that in DBI inflation, the scalar potential dominates the energy density even for relativistic 
rolling. It can then be seen that the term on the right hand side is a decreasing function as one 
approaches the relativistic limit. This reduces the solution space for A^ so that it approaches 
zero as 7 — )■ 00. 

Finally upon variation of the action we find the coupled Maxwell equation 

d, [V^Fr{<l^)Fl{<t>)F>^^) = 0, (3.12) 
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which mixes the infiaton with the gauge field. Inserting our ansatz from section 2.1 we find that 
the equation becomes 

■ / F{ ■ 2Fl ■ Fl<j)^ 2F2U \ 

^, + ^,^„ + 4. + ^, + ^* + ^j-^ + ^^j=0. (3.13) 

Solutions to the Maxwell equation include the following - similar to that in the Bl-theory 



PAe 



-a—Aa 



Ax ^ o— , (3.14) 

where we denote the constant of integration as Pa- In this case Pa is a measure of the charge 
carrier density on the world-volume because the gauge field arises from excited states of open 
F-strings that end on the brane. 

In the non-relativistic limit (7 ~' 1), the scalar potential will effectively dominate the right 
hand side of the constraint equation (3.11). Therefore for a sufficiently 'large' potential, this con- 
dition can easily be satisfied. In the ultra-relativistic limit, the right hand side is proportional to 
7~'^(y — 7F1/2) which may not be very large without fine-tuning, thereby making the constraint 
equation (3.11) very difficult to satisfy. Clearly the gauge field will have the most dramatic effect 
on relativistic inflation, but may also play an important role in the non-relativistic limit. We 
can again, like in subsection 2, adapt the constraint on the energy density during inflation, in 
the form of a ratio of energy densities to be approximately constant (cf., e.g., (2.17)). In the 
case of the DBI theory, the new expression for the ratio takes the form 

^=2FiF|(Fi(7-c?) + y)' ^^-'^^ 

which clearly relates the various parameters in the model. Examining the ratio ()3.15p for suffi- 
ciently small values of - and considering the critical limit where the gauge field energy density 
is almost constant during inflation - we find the constraint 

FiF| ~ e-^^ (3.16) 

which is the generalised extension of the known results from Einstein-Maxwell theory. To further 
assist in understanding the role of the gauge field, we use the parameter c (cf. Appendix A), 
which conveys the deviation from the above result. More precisely, the solution we must consider 
is FiF^ ~ g-4ac _ -^i^igj^ ig i]xQ most general expression. For c > 1, the gauge field terms will 
rapidly come to dominate the dynamics at late times. 

Before proceeding and specializing with different settings of our DBI framework, let us 
indicate that the anisotropics are governed (at leading order) by equations (j3.7p . ()3.8p above. 
Under the assumption that they are small, and obey slow-roll behaviour, we can immediately 
write down the relevant Hamilton-Jacobi equation which describes their dynamics 

da^\FiFiJ\ 67 JF^{J-q) + V■ ^^ > 

In the non-relativistic limit, the 7 terms reduce to a constant (given by 1/3), and the magnitude 
is then set by whichever term dominates the scalar energy density. For ultra relativistic motion, 
the middle term appears to be linear in 7, and if the scalar potential is dominant this suggests 
that the anisotropics may be large. However if the F17 term dominates the energy density, then 
the 7 factors actually cancel - and the anisotropics are then determined by the new term F^F^. 

14 



3.1 Non-Relativistic limit 



The non-relativistic limit of the theory emerges when 7 ~ 1 + 2-^2</'^ + • • •• Although this 
is a modified version of slow-roll inflation, there is a non trivial coupling to the gauge field and 
therefore one may expect back-reaction to be an issue. Note that for q = 1 the effective potential 
Veff = V{(j)) + Fi ((/))(! — q) reduces to the scalar potential due to supersymmetry, and when 
q = we have a purely non-BPS system - which also allows for V = and the dynamics are 
dictated by Fi alone. 

Expanding the DBI Lagrangian to leading order yields a corresponding effective field equation 

+ 3a,^ + 0-- — --- = 0, (3.18) 

where herein dots denote time derivatives, primes are derivatives with respect to the scalar field 
and we employ the following variable definitions 

Q ^ FiF2 (l + ^F^^Ff^A = F1F2 (1 - I) , (3.19) 

W = Fi(g-l)-y-Fi|, (3.20) 

where x = — ^F^yF^'^ . We have not made any assumptions about the background geometry 
at this stage, and this result is therefore quite robust. Let us study the background evolution 
by starting with the isotropic inflationary limit and initially set the gauge field and (p terms to 
zero, in which case the above equation of motion reduces to 

■ d? d (q-l) d V 

which is reminiscent of a canonically coupled scalar field solution - except that there is an 
additional driving term proportional to (jP'. Note that BPS configurations (g = 1 in our language) 
simplify the equation of motion significantly since the second term in (3.21) vanishes in this limit 
(with the usual assumptions about regularity of the background function) . The troublesome (fP' 
terms, which have no analog in canonical slow roll models, will always contribute - unless we can 
consider limits where F1F2 is constant. This amounts to localising the solution on a curve in 
parameter space, and is what we consider in this paper. General solutions to the above equation 
will be explored in the future. 



3.1.1 F1F2 Constant 

With the assumption that F1F2 is constant, the field equation becomes a modified version of the 
canonical scalar field equation. We can immediately write down the scale factor as a function 
of the inflaton using the Hamilton-Jacobi formalism, for an arbitrary scalar potential 
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For a simple chaotic inflationary potential such as ^''^'^4''^ theory, we can then see that the 
back-reaction constraint (|3.16|) implies the following general relationship 

F,F, ~ exp [-^^) ■ (3.23) 

which dictates the functional form of -F3, given that Fi is fixed. 

We now proceed to include the gauge field terms in the modified Klein-Gordon equation 
which we write as a perturbative expansion in x, and takes the general form 

■ / 2x\ 6^ d , , V / X 

(x-(x + l)(g-l)) a , ,^, , x(F2<A^-2) a , ,^^, ._.. 

+ 2F, 9^^"(^^^+ 2F, a^^"^^^^^)' ('•'') 

although the (p'^ term in the last bracket in the expression above is actually subleading in the 
non-relativistic expansion - and can be neglected. Note that the source term also appears to be 
coupled to the scalar potential, which is due to the non-canonical nature of the theory. To solve 
this equation, we use the energy density constraint (j3.15p to solve for F^ under the assumption 
of a quadratic scalar potential - since this is the simplest analytic solution. One can solve this 
equation numerically, however we present only analytic solutions in this paper - leaving a more 
exhaustive analysis to future work: 



(a) Let us initially consider Fi to be constant, which (by assumption) implies that F2 is also 
constant. Using the back-reaction condition to solve for F^ we can then find the solution 
to the scalar field equation in the presence of the gauge field. At late times the solution 
can be calculated to converge to the following 

cPJ(6F2F2C-m2Mp2)' ^ ^ 

where the right hand side is clearly a constant which vanishes for c = 1, therefore the 
electromagnetic energy density will also be constant. We see that FiF^ ~ l/a^ in the 
regime where Pa is constant, but more importantly it highlights the fact that c ~ 1 as the 
back-reaction becomes relevant. 

For a canonically coupled field, there is attractor behavior since the gauge field energy 
density exhibits tracking behaviour [18]. In our case, due to the non-canonical nature of 
the action, the source term in the field equations also couples to the scalar potential. If we 
demand that the scalar and gauge fields are effectively the same magnitude, so that they 
are both source terms in the field equation, then we find 

2,^2 xc{FiF2)^ / , X 



at leading order. Inserting this into the relation (jS.lSp we discover that 

7^~— r--f— -. (3.27) 
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Given that inflation occurs for ~ O{10)Mp in tliis model, we find 7^ ~ xlO~^ /{F1F2), 
which can be vanishingly small for sufficiently large F1F2. This ensures that the gauge 
field energy density is negligible when compared to the scalar energy density and therefore 
we can consistently neglect its contribution to the Hubble parameter. Note that this is 
almost the same result as obtained in canonical models, aside from the factor of F1F2. 
This ensures that 7?. ~ xlO~'^ /{F1F2) is an attractor solution. If 7^ is initially smaller than 
this, the behaviour of the scalar field drives pA to increase (with c > 1), reaching this value 
from below. For configurations where this quantity is larger than the attractor solution, 
the inflaton can climb back up the scalar potential (due to the source term in the equation 
of motion). Thus pA is forced to decrease rapidly, and the attractor is reached from above. 
This confirms that even this non-canonical model we expect tracker behaviour. 

As far as the anisotropies are concerned, dropping the a term in the anisotropic equation 
of motion (|3.8p , combined with the back-reacted solution implies the following 

a 2Pl 6(c -l)Fi m^M^ 

^ " 3m2(/)2 cPl{i6FlF2C - m?M^) ' ^'^''^^' 

indicating that anisotropies increase during this accelerated phase for c > 1 and QcFfF2 > 
m?Mp. Of course, in order to determine this result we assume a perturbative a as before. 
Note that, as in the case of Born-Infeld theory, the increase in anisotropy is determined 
by the scalar potential. Both solutions (j2.14p and (j3.28p increase like l/4>'^ - although in 
the latter case the anisotropies also vanish when c = 1. 

(b) Consider now a regime of solution space where the logarithmic terms dominate in the back- 

reacted equation of motion - still assuming that F1F2 is constant. The solution for the 
scalar field in this instance becomes 

„-4a 2 7\//'2 

g_4a-FiF2C<^2/A/| ^ "" ^^> (3 29) 

nF2F|c2p2(l _ e-4a) + ^2Af2gFiF2C0g/M2 ' 

where n is an integer depending on whether Fi{n = 1) or i<3(n = 2) is fixed to be constant. 
Note that this is similar to the result obtained in (|3.25p with the additional dependence on 
the scale factor appearing on the right hand side which ensures the solution decreases as a 
function of time. The scale factor dependence arises precisely because of this logarithmic 
running. The solutions are valid in a regime where the scalar mass in the quadratic 
potential satisfies the following bound 

9 nxcF?F2 / x\ , , 

™^«-jjf^(l--) (3.30) 

where n is defined as above, and when n = 2 we must recall that Fi depends on the 
inflaton. This bound can be satisfied for a (small) region of solution space and is therefore 
a physical solution. If the parameters are chosen so as to satisfy this bound, then one 
can easily show that for initially small anisotropies, these are actually decreasing during 
the inflationary expansion - indicating that the isotropic universe would be a late time 
attractor. 

(c) It is also of interest to explore a solution branch where we fix the functional form of -F3, 

but with the same assumption that the logarithmic terms are dominant with respect to 
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the scalar potential. This should then modify the above result for the scalar field. The 
ansatz we select is F3 ~ Lcj)"" - which allows us to find the analytic solution 

which is in fact identical to the solution above for n = 1. This suggests that the anisotropies 
are decreasing in this regime - since their overall magnitude is a decreasing function of a. 
Let us therefore consider a solution where F\ ~ Lcf)'^ - which continues to fix F2, but now 
with F3 as a constant. The resulting expression for the back-reacted field equation is 

^'(«) ^ ^£l£^g-4a-FiF.c<^V2Af,^ (332) 

which admits a complicated solution of the form 

2 



-cF^F^cjy /Mp _ g^p 2InverseErf iErfi 



Mp 



F1F2C 



^(FlF2c)3/2p2(l_e-4a 
m^Mp^/27r 




where Erfi is the imaginary error function [3T] and therefore one sees that even for c = 1 
there are non-trivial contributions to the anisotropy equation. 

(d) We have considered the equation of motion where the logarithm terms were dominant, but 
more physically for the non-relativistic case, we can anticipate the scalar potential to be 
the largest contributor to energy density. We may again choose to set either Fi or F3 fixed 
to be constant, which leads to the resulting expressions shown below; 



-Aa-cFiF2<f?/Ml 2Fi(l - c) 



{F3 
A 
(1 - 2c)e-2" 



PaFsV^ 



(Fi). (3.33) 



The first solution is constant for constant F3 and clearly vanishes for c = 1. The second 
solution is for constant Fi, but clearly decreases with time. Moreover the solution does 
not vanish for c = 1, rather it vanishes for c = 1/2 - indicating that anisotropies will be 
increasing in this case. 



3.1.2 AdS^ Solution 

As an example of a particular solution, let us consider the case of a pure AdS^ embedding, 
generated by N coincident i53-branes at large A^. The background functions must now satisfy 
the following expressions Fi = cp^/X^jF^ = l/y/Fl and F2 is constant. Since the background 
functions are explicitly known, the gauge coupling is fixed - in fact it is unity which leads to a 
trivial result. The solutions depend explicitly upon the particular brane being embedded into 
the theory. 



18 



We focus initially on non-BPS configurations where q ^ 1- In this case, we must ensure that 
A >> 1 for analytic solutions. We will again assume that the scalar potential is quadratic. With 
the gauge field set to zero, we find the equation of motion at leading order in A becomes 



4>'{a) 



2Mp2A2 



+ ..., 



(3.34) 



where the higher order terms are sub- leading in A^. 
solution 



The above equation admits the following 



/(«) ~ </>; 



^2 



(3.35) 



where we neglect the sub- leading terms in A. This leads to a cancellation of the scalar field mass 
in the equation of motion, therefore the solution does not appear to depend on it at leading 
order. The corrections coming from the gauge field lead to the following master equation 



</.'(«) 



2M^X^ 
' F205 



1 



A^P^e-^" 



(3.36) 



which must be solved perturbatively. Let us define (p = x + Pa^ + • • • at leading order - where 
the solution for x is the one found above in (j3.35p . The back-reacted solution therefore has the 
leading order solution 



/ ;^2g-F20[>/(3A2Af2) 

(«)~X 1 + ^, <E 



'3A2M2 



4e 



2 
3' 



Fo 



3XUII 



(3.37) 



where we use the short-handed notation E[x,y] = ExpIntegralE[x, y] for simplicity |31| . Note 
that the boundary conditions ensure that (p ^ (f)Q at a = 0. Since both Fi and F^ are defined 
by the background geometry in this instance, we may find it hard to satisfy our gauge density 
constraint (|3.15|) . Indeed, one can see that 7^ only depends on the ratio e~^°/(/)2 which cannot 
be constant at this order of approximation. 



In the case of g = 1, which is the BPS configuration (in the static limit) we now need to 
assess whether we keep the (p"^ terms in the equation of motion. If we set them to zero, then 
we recover the expression (j3.34|) as an exact result, not a perturbative one. The back-reacted 
solution will also then follow trivially. Instead let us keep the quadratic terms in the equation of 
motion. We can then solve this equation in the absence of a gauge field to obtain the expression 



3a4> 



+ 



(3.38) 



where we must ensure that the following condition is satisfied 



1 » 



3F2(/)6 
16A2M2 ■ 



(3.39) 



This can be done if we also assume A^ >> 1, which is good for the supergravity approximation. 
In this limit of the theory we find that dynamic solutions are exponentially decaying 



/)oe 



-3a/2 



(3.40) 



19 



and if one includes the backreaction at leading order, the solution is 

Ha)^e-''^/'{cPl±X'Pl{l-e'nf\ (3.41) 

where we must again impose the constraint condition ()3.39p . Intuitively these results make sense, 
because we are forced to consider a perturbative expansion in the mass term, which effectively 
decouples the scalar potential from the theory. The scalar field is driven purely by the potential 
generated by the AdS background. 

The energy density ratios in the two cases of interest may be written 

where TZi is the solution arising when we neglect (j)^ terms in the equation of motion. Note that 
TZi is initially very small, but increasing during inflation indicating that the energy density of the 
gauge field is highly suppressed. As (j) continues to decrease, corresponding to the limit where the 
probe branes are nearing the D3-branes, the ratio rapidly starts to diverge indicating that the 
gauge field contribution overwhelms the scalar energy density and back-reaction dominates. The 
solution for TZ2 also increases as the brane approaches the stack, although one must be careful 
to ensure that the parameters satisfy the constraint ()3.39p . One notes that, when compared to 
TZi with similar choices of parameters, that 7^2 is significantly smaller in magnitude than 1Z\, 
and the gauge field domination occurs at later times. 

Regarding the anisotropics we see that herein the equation of motion in the BPS case g = 1 
(dropping the (fP' terms) can be written as follows; 

^ ' ^" ' , (3.43) 



r\j 



which are increasing rapidly in this instance due to the exponentially decaying behaviour of the 
scalar field. In the case where we neglect the cjP' terms, the solution becomes 

. 4P>V,_£M^V (3.44, 



n^ 



a 3mV^ V 3A2M2 

at leading order in A for all q. The result is that the anisotropics are increasing during infla- 
tion, recalling that the parameter space is tightly constrained, and at a smaller level than the 
previously considered solutions. 



3.1.3 Tachyonic solution 

Let us consider the non-relativistic expansion for the tachyonic solution in a non-BPS conflgura- 
tion, where F2 is a constant, Fi is the tachyon potential and g = 0. We will take the potential to 
be of the form Fi ~ Vo/cosh((/>/L), where L is an unknown dimensionful parameter. We leave 
this arbitrary because it is highly likely that the tachyon potential in curved space is different 
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from the one derived using BSFT (Boundary String Field Theory) in flat space [j- Additionally 
we will also set the scalar potential to zero, ensuring that the inflaton dynamics is driven only 
by open string condensation on the world- volume. The resulting ratio constraint from ()3.15p 
becomes 

F^Fi ~ Z^e-^""", (3.45) 

where Z is a dimensionful parameter of the theory (cf. Appendix). 

The background field equation with the above tachyonic potential results in the following 
expression for the inflaton 

.a„h(|-)~exp(5^). (346) 

where we have dropped the constants of integration. It is well known that tachyonic inflation 
requires severe fine-tuning, however the literature has have always assumed the validity of the 
flat space potential in curved space. Indeed when warping is taken into consideration, this fine 
tuning decreases allowing for potential inflationary trajectories |50j . In our case the warping is 
essentially embedded in the definition of L, which can lead to inflation when it is sufficiently 
large since it suppresses the mass of the inflaton. 

The gauge coupling in this instance can be written as follows 

FiFi ~ T^e-^"'^ cosh (2arctanh (e"'^)) , (3.47) 

where we defined r = M^/{2F2L'^Vq). Note that this coupling is non-positive, and increasing 
towards zero from below. This is true even if the brane was a ghost brane, with negative tension. 
The backreacted equation of motion can be written as 

M^ / 4.t\ f) 
which admits the following asymptotic solution 

*"^^^ UJ ~ MMi^ "' [ 6ic-mL^voz^ ) ' ^'-''^ 

where we have included the initial boundary condition on (p such that (f> ^ e « L at a = 0. 
The gauge coupling is, in general, a rather complicated function of the scale factor 

F,Fi ~ ^^ cosh (2arctanh (tanh (^) exp (^-^^j^—^ (Pi(-1 + e^^^^^^^) + 3aZ2(c - 1)))))) 

(3.50) 
for arbitrary values of c; However, for c = 1 the solution reduces to 

r-,2 —4a 

FiFi ^ cosh (2arctanh (tanh (^) exp (^[4?! + SZ^]))) , (3.51) 



^See the recent paper [49] for additional clarification 
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where r is the same quantity defined earher. Because of the perturbative correction in P^ the 
gauge couphng for c = 1 can be a positive, decreasing function (as a function of the expansion). 
However due to the algebraic structure of the solution, there is a critical value of the scale factor 



3Z 

a, 



Infcothf^)) (3.52) 



r(4p2+3^2) V \2LJ 

at which the solution is singular, therefore we can only trust the solution for a < Uc- When the 
background is uniquely specified in this case, we see that decreasing solutions require r << 1. 
When we initially fix r, we then see that a decreasing gauge function requires Z'^ << Vq. 



3.2 Relativistic limit 



The relativistic limit occurs when 7 >> 1, corresponding to <j)'^ ~ ^/F2-, and the scalar 
equation of motion takes the following form 

.. F-M (1 - f ) (3c. + ^1 ln(F.) + ^'^ + ]'^^^'' I ln(F.) j (3^53) 

+ ^.7*. (2d + 4 HF.F,)) + i (l - f ) I ln(FO - ^| HF^F,) + 11^^. 

Since 7 is controlled by F2, and we know the relation between this function and the inflaton, we 
will treat this as an unknown variable allowing us to eliminate the velocity from the problem. 
The equation of motion can then be solved once we specify the potential and the function F\. 
Unlike the non-relativistic limit, there is no 'trivial' simplification that one can consider. The 
best approach turns out to be fixing the potential to be quadratic, as before, and assuming that 
Fi satisfies power law behaviour, Fi ~ Lcjf where L is herein a constant with mass dimension 
(4-p). 

In general, our assumption that 7 >> 1 allows us to solve the system explicitly once we 
specify F2. To illustrate, let us assume that F2 ~ W(j)^ . We can then integrate the equation of 
motion directly to obtain 

-^(^^ ~ ^ /fi/ , 4^ (^^tf '"'^^' ('Ao -ct>)+ m2 + p)) , p^-A (3.54) 



Mp y/6{p + 4) 
where W, T are constants. We can then investigate the following cases: 



p = -A (3.55) 



(a) Let us assume that V >> Fiq as a constraint on the Hubble equation. This is the limit 
in which (canonical) DBI-inflation can occur [5H I52j. Physically this corresponds to the 
inflaton being strongly damped by the Hubble factor, so that even though it is relativistic 
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in velocity, it does not travel a large distance in field space. Knowing the scalar field 
solution, we can then use the Hubble equation to determine the scale factor 

«W ~ ^^-j^^^^{VT^^t'^/\^o-<P) + tcP{2+p)), p^-4. (3.56) 

"" '^^ ' ~^" ' p=-4 (3.57) 

where we have assumed that (j) is non-negative for simplicity - and imposed boundary 
conditions on the scale factor so that it vanishes at the start of inflation. 

One could include the gauge field corrections to the above results, by employing the relevant 
term in the Hubble parameter. However it will turn out to be much simpler to use the 
scalar field equation to solve for F2, and then compute the back-reaction on this variable. 
We will now consider this strategy where F2 is unknown and Fi has power law behaviour 
Fi ~ L(j>P. The general equation is quite difficult to solve analytically once one includes the 
back-reaction. Therefore we illustrate the results for the case p = 2, which has a leading 
order expansion of the form 

^'" 3m204 (^ 3L02(473(c-l)(zy2 + 2)-3i/2) 

(3.58) 
where the ± sign arises from the particular choice of sign in the Hubble equation, and 
z/2 = -^. One notes that the leading order dependence on the inflaton is l/</>^ and 
therefore one expects that the scale factor should initially run logarithmically (at leading 
order) 

which can rapidly become large due to the overall pre- factor of 7^, and the constraint that 
V » Fi7 always ensures that V2 » 2. 

We can immediately ask what happens to the anisotropies in this limit. Indeed, in the non- 
relativistic limit we discovered that they increased during the inflationary epoch because 
of the back-reaction. A quick calculation in the relativistic limit implies 



2 /J, \47^(i'2-2)(c-l)/3 



where c is the parameter arising from the gauge field condition FiF^ ~ e~'^°"^. One can 
therefore see that anisotropies increase quite rapidly in this instance - significantly faster 
than in the non-relativistic limit due to the dependence of 7 in the exponent. 

(b) The converse limit, where -F17 >> V is also interesting, since - like the tachyonic theory in 
the non-relativistic limit - the gauge field constraint in (j3.15p imposes the condition that 
FfF^ ~ Z'^e~^^'^. We will again assume that Fi is power law to simplify the equations of 
motion. The isotropic theory admits the following solutions 

Ml{^\2 + p)-2pf 
^^ Vl^^fL^v ' ('•'') 

a ^ (7!(2±z)^Mi,fM, (3.62) 

6 V / 
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which are significantly diflFerent from the potential dominated regime. In particular we see 
that the scale factor is logarithmic for all values of p. This solution actually has many 
similarities to the non-relativistic tachyon theory, since we can again write the gauge 
coupling function as follows 

Including the gauge field back-reaction, we then obtain the leading order solution 
^^ ^ Ml{^^{2+p)-2pf / ^ .^y(c-i)(.^(2+rt-2p)/3 (8c + 3±l) ^ \ 



127L02+P I Z2 V (/) 7 (472(0 - 1) =F 3) 

(3.64) 
where the correction term will clearly be constant when c = 1. The choice of sign again 
arises from the Hubble expression. Upon integration we then find the following solution 
for the scalar field 

^T l^ + 4Z2(c-l)AA472(c-l)^3j' ^^-^^^ 

which simplifies considerably when c = 1 to become 

where we have herein used the notation N = 72(2 +p) — 2p. This latter expression can be 
inverted to obtain the scale factor in terms of the Lambert W function [31], which admits 
the following expansion; 

NZ2 / 6 \ „ 12±1 



at leading order in P^. The solution for the scale factor therefore appears to be non- 
perturbatively corrected by the presence of the gauge field. The gauge coupling then takes 
the schematic form 

FiF| ~ a-Pe-^"^ (3.68) 

which exhibits a minimum if Oc = — p/(4c) where p,c have opposite sign and c < - 
however since we require c > for physical solutions, this is most likely not a physical 
result. However this point is a maximum when p < with c > 0, thus there are solutions 
where the gauge coupling initially grows in strength, before reaching its maximal value at 
which point it starts to decrease. This unusual behaviour arises from the analytic structure 
of the DBI action, and therefore has no natural analog in terms of canonical scalar field 
models. 

(c) Let us briefiy discuss the tachyonic solution in the relativistic limit. The isotropic equation 
of motion yields the following solution for the scale factor valid at large 7 

a = -ln(cosh(^ j j . (3.69) 
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Including the leading order back-reaction we then find the following solution 

(3.70) 



cosh ( y ' 




p2 -,l/A\mc-l) 

f_A I2a{c-1) 



in terms of the Lambert W function, and valid for all c 7^ 1. This modifies the gauge 
coupling function resulting in the expression 




p2 nl/4Xl/(-l) 



Z2 



Agl2Q(c-l) 



(3.71) 



which is not necessarily a decreasing function. For c < 3/4 we see that this function 
actually increases with scale factor, whilst for c > 3/4 the function is monotonically 
decreasing, with an amplitude set by the ratio Z^/Vq. For sufficiently small electromagnetic 
energy density we can expand the Lambert function and we then find that 

F,Fi^^{-7nHf+^, (3.72) 

which corresponds to a spectral index of n = 7 + c for the density ratio, which should 
be contrasted with the tree-level result which yields n = 1 -|- 4c. Therefore we see that 
even a tiny electromagnetic field in the tachyon case leads to additional dependence on 
the expansion parameters, and a significantly larger magnetic field for c < 2. 



4 Discussion and Outlook 



In this paper we have investigated a class of corrections to inflationary solutions arising from 
the introduction of an electromagnetic field in a non- linear context. We first considered the pure 
Einstein-Born-Infeld theory, and found that, in order to have satisfactory inflation, we needed 
to promote the BI coupling to have scalar field dependence. For power law solutions we then 
used the WMAP data set to restrict the functional form of the coupling, subsequently retrieving 
a range of parameters inducing models observationally consistent. 

Following this, we employed a generalised approach to D-brane inflation where space-time 
backgrounds |j were not necessarily asymptotically AdS (cf. [25] for an analysis including an 
asymptotically AdS setting.) This allowed for more richer classes of couplings than in pure 
Einstein-Born-Infeld settings. In all cases, the presence of the gauge field was subleading in the 
corresponding inflationary trajectory. Moreover, anisotropics tend generically to increase with 
time; The anisotropics were treated perturbatively in this paper, but a more detailed analysis 
in such models would be welcome. The one case where this was not occurring was when the 
logarithmic terms dominate the non-relativistic equation of motion: The anisotropics decrease 
with time, provided that the scalar potential satisfies a (stringent) non-trivial bound. 



^Such backgrounds are generated by gluing warped throats onto a Calabi-Yau three- fold. 



25 



Concerning the analysis above summarized, it is of relevance point out that the scalar poten- 
tial was assumed to be of a simple power law form. Whilst this is convenient for comparing so- 
lutions to the canonical (non-relativistic) limit, it is not necessarily easy to obtain such solutions 
within a string theory context: Although quadratic potentials are common within supergravity 
theories, stringy instantons typically give rise to exponential potentials; Therefore some of the 
results obtained in this paper will not be valid in a more complete string theory embedding. 

Finally, regarding the evolution of magnetic seed fields in a DBI context |j, a brief analysis was 
provided in Appendix A. The essential feature is the new couplings that the DBI configuration 
induces. It was seen that the electromagnetic coupling depends on two parameters. We chose 
to fix one of these parameters by demanding that the energy density ratio was constant. Two 
classes of solution could then be identified. The first was a generalised version of the canonical 
coupling, given by e~^"'^, which led to scale invariance when c = 3/2. The second class of 
solutions were specific to the model under consideration. Notably the tachyon (relativistic and 
non-relativistic limits) and the ultra-relativistic limit of the D3-brane theory. In such models the 
scale invariance arose for values of c < 1, as demonstrated explicitly in the case of the tachyon 
which has a coupling of the form e°'^^~^'^' . For special classes of DBI models in the relativistic 
limit we found that the coupling behaved like a~^e~'^"'^ where p arose from the assumption of 
power law dependence. Overall, we found that the obtained magnetic field spectrum is typically 
larger than the observed bound unless it was created primordially (with subsequent amplification 
via the dynamo mechanism). 

As a last note, allow us to indicate that a maximal bound on the strength of the gauge field 
could be generically established. When saturated, implying that a non-zero field is generated at 
sufficiently small scale, a residual cosmological constant emerges - independent of the form of 
the scalar potential. It is not, therefore, inconceivable that a theory could eventually contribute 
to resolve the the dark energy problem (and possibly the issue of primordial magnetogenesis) . 



A Magnetic field generation 



We have explicitly considered a non-zero electric field component as a solution to the field 
equations (|3.12p - (j3.14p in the main body of the paper. For the DBI theory of L)3-branes this 
corresponds to turning on F-string flux on the world-volume, and is therefore rather natural. 
The magnetic field in this case corresponds to the Dl-string flux. The Z?-brane theory exhibits S- 



®The magnetic field on the D3-brane can be seen to be generated by exciting the open string degrees of freedom. 
Since the F and D strings are S-dual, the most general configuration would be a mixture of both string solutions. 
Since the increase in world volume flux tends to increase the mass of the moving brane, it is likely that the brane 
will 'sink' lower in the throat therefore making inflation harder to occur. However there may well be a set of 
backrounds where this does not occur. Turning on such flux, however, could be useful for transferring inflationary 
energy to the standard model sector (should they be localised on seperate branes), because the open strings on 
the inflationary D3-brane could attach themselves to the SM branes at the end of inflation. Such a process would 
correspond to the direct transmutation of energy from inflaton sector to the standard model degrees of freedom. 
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duality invariance, which encompasses the electro-magnetic duality of ordinary Maxwell theory. 
Consequently, one can deduce the levels of the magnetic field from the electric field, and vice- 
versa. 

Let us therefore consider the manifestation of such fields during inflation. In general this 
is a highly technical problem, hence we must consider simplifying the background to make the 
analysis tractable. This is usually done by initially consideration of the flat FRW geometry (cf. 
section 1), and moving to a conformal gauge such that t = f e~°'dt. The gauge field can then be 
expanded into annihilation/creation operators and mode functions A(r, k) satisfying the Fourier 
space equation 

^ + ^ln(FiF|7)^ + :^A = 0, (A.l) 

where k is the comoving wavenumber. If we now introduce another change of variables rj = 
J dT^~^, then the above equation becomes the more familiar wave equation 

-^ + 2- HF.Fl)-^^ + k^M,, k) = (A.2) 

and we see that 77 = —e~°'/{'^H) for constant ^H. This definition allows us to divide the mode 
functions into sub and super horizon modes, where the former occurs for k » ^He" and the 
latter for modes satisfying k « ^He'^, for some characteristic time scale ry^ ~ —1/k. Making 
the following identification 

Y' _l dHF,Fi) 

Y = 2 d^ ' ^^■'^> 

using primes herein to denote derivatives with respect to rj, we obtain the following Schrodinger 
type expression through the change of variables v{rj, k) = Y{rj)A{ri, k): 

= v" + v(k^-^p\ (A.4) 

Sub-Hubble modes admit a solution given by the WKB approximation, written in terms of the 
gauge field mode expansion A 



g-ikrj 



111 / TT 7 

y/2kF\Fi 
whilst the super-Hubble modes admit a solution of the form; 

drj 

F^Fr 



(A.6) 



out - Ci{k) + C2{k) —^, (A.7) 



where the Ci are constants of integration which can be determined by matching the in and out 
modes (and their derivatives) at horizon crossing He^ = k. Note the above expression is valid 
at leading order in k. Neglecting the subsequent decay mode for in Aout one can then write the 
following general solution for the mode expansion - following the arguments presented in [25] 



|A(r?,A:)| 



2 1 



r\j 



2k{FiF, 



2\ 






2k{F,Fi)J L (FiF|) 



(A.^ 
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where the asterisk subscript denotes that the quantity is evaluated at horizon crossing. Assuming 
that 7, H are both constant during inflation, which is a good approximation, then we can use 
the foflowing identity 

d, = -^ (A.9) 

to simphfy the above expression. One can calculate the vector field correlation function to obtain 
the electric field power spectrum, and using the inverse Maxwell relation Bi = eijkF^^/2 we then 
see that we can write the (proper) magnetic field power spectra in the following manner 

|5|2~2/fc2e-4"|A|2, (A.IO) 

which has the correct behaviour in the radiation dominated epoch so that B ~ e~^". The 
derivation of this expression relies on the fact that we assume instantaneous reheating, and 
that the conductivity becomes much larger than H immediately after infiation. Introducing the 
density parameter Q{r],k) = Pb/p^ where p^ is the energy density of radiation, and ps is the 
gauge field energy density (per unit logarithm), we have the definition; 



0(r?,A:) 



15 (ke-'^^Y {FiF^) 



27T^N V Tr J (FiF^) 



1 _ ( , + (^1^3 )'* A . r iFiFh*H 



2k{FiFi)J ^ (i^K 



2 



(A.ll) 



where A^ is the number of massless degrees of freedom at reheating, and Tr is the reheating 
temperature (a subscript R denotes the quantity evaluated at the reheating epoch). Assuming 
that the gauge coupling exhibits power law behaviour such that 

FiF| ~ const X ( — J , (A.12) 



we can then see that 



indicating that the spectral index is given by 4— p. The precise value of p can then be determined 
once the precise theory is specified. 

At scales of order L = 2t: /k [Mpc], it was subsequently derived in [2S] the following bound 
for the present magnetic field assuming a reheat temperature of the order IQ^^GeV, N = 100 



where we have defined the gauge function as FiF^ = F{r]) and W{p) is a p-dependent function 
that diverges as p — )■ —1 and asymptotes to W ^ 0.25 as p — )• ±oo. Note that we must have 
\B\o > W~^GeV to be consistent with observed field strengths up to 1-Mpc without envoking a 
dynamo mechanism. Clearly this is a very tight bound for all models to satisfy. The expression 
for the field at decoupling (a^ec ~ — 31n(10)) is essentially the same. The only difference is that 
we have a shifted overall exponent 10^^^^^^ and must of course replace F{0) by ^(dec), with 



28 



which must now satisfy the bound |i?|dcc ^ 10~^'^G for cosmic seed fields. This is a less stringent 
bound because in this instance we can assume there is a dynamo mechanism which serves to 
amplify the field at late times. 

In terms of the gauge field coupling function -F1-F3 , we can determine two immediate situa- 
tions of interest (cf. eq. p.lSp ): 

• Firstly, there is the case FiF^ ~ e"^"'^ for both relativistic and non-relativistic cases, 
where c is a constant. This occurs when the scalar potential dominates the total energy 
density of the model. This is similar to the condition established in the BI section and it 
can be shown that the mode function after horizon crossing is proportional to ei'^'^-^)'^ - 
indicating that we must consider positive c, because negative c leads to a decreasing gauge 
field component, which we naturally expect to have a smaller contribution to the total 
energy density. Inserting this expression into the mode expansion of the gauge field, and 
obtaining the result for the proper magnetic field, one sees that (at the end of inflation) 

|i?|2~/?fe2"/(2c-3)^ (A.16) 

where Hj is the value of the Hubble parameter at horizon crossing. For a flat spectrum we 
therefore recover the constraint that c = 3/2. Indeed one can show that this result implies 
that the energy density of the electromagnetic sector rapidly becomes comparable to that 
of the inflaton. For the marginal case where c = 1, we can explicitly solve the gauge field 
equation of motion. Matching the solution to the vacuum in the sub-horizon limit, we can 
identify the integration constants, and write the solution in terms of Hankel functions 



/ „(l-4c)/2 , ^ /Ar - 1 \ 

which is plotted in Figure 2. In terms of magnetic field spectral index and written as 
n = 4 — p, we see that this immediately implies n = 4(1 + c). Since this is positive definite 
for c > 1 we find that the spectral index is highly blue-tilted and damped for larger values 
of c. Note that c = 3/2 yields an approximately flat spectrum as anticipated. 

The second case of interest occurs when Fi{'y — q) >> V, which could arise in the ultra- 
relativistic regime where Fi'j is very large. The constraint equation (|3.15p then implies that 
F^F^ ~ Z'^e~'^""^^, where ci is a dimensionless constant and Z is a dimensionful constant. 
Since the gauge coupling depends only on FiF^, we must specify one of these two functions 
completely, before being able to determine the magnetic fleld spectrum. This new limit 
is only possible due to the non-linear nature of the DBI action, but appears to depend 
explicitly on one of the background functions - unlike the previous case. Interestingly this 
is precisely the combination of parameters that determines the scale of the anisotropics in 
(|3.17p . Our emphasis has been on determining F3 for a given value of Fi unless otherwise 
stated. Furthermore: 

— Empolying ()3.63p induces solutions of the form 

A{k,v) ~ ^/^"«^)/' (^iJ (^^'^^) +^2Y (^l^,fcr?)) , (A.18) 
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Figure 2: Figure showing Re{A)k^^~'^'^'''^ as a function of z{= —kr]) for different c. We have also 
absorbed a factor of {'^H)'^'^ into the gauge field for simplicity. The level of damping increases 
with c. 



where Y and J are Bessel functions j31j . The assumption of large 7 ensures that Qi 
tends to 4c from above for all values of p > 0. In order to obtain smaller values of Qi 
we must ensure that p is negative - in particular one can see that the 7 dependence 
drops out for p = —2 which is plotted in Figure 3, although note that Qi ~^ ~4 
asymptotically. The solution simplifies in the case of Qi = 1, which occurs when p 
takes the critical value pc, 

which asymptotes to pc = —2 for large enough 7 independent of the value of c. 
Therefore, if Qi is constrained to be unity, then we see that p — )• —2 is an attractor 
solution. Normalising the solution using the sub-Horizon modes, we can extract the 
coefficients arising from integration and write the solution: 



A. ~ e-^^'3i/4 



7? 



i-Qi 



2(-7F)Oi 



^(2) 



Qi 



1 



-,A:ry 



(A.20) 



where H^'^> is the Hankel function of the second kind, which is illustrated in Figure 
(3). Note that the attractor solution p = —2 ensures that the gauge field is strongly 
damped. Other values of p tend to lie on the same curves, indicating that the solution 
is insensitive to the precise form of the power law. 

From ()3.69p . we can use it to solve the gauge field energy density constraint to obtain 
the electromagnetic field as before. The solution is again a Hankel function, but if 
we solve for the magnetic field we obtain the following scaling 



which indicates that c = 3/4 yields a fiat spectrum at the end of inflation. 



(A.21) 
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Figure 3: The real part of the expression Avk^ ^ plotted as a function of z{= —kr]) for c = 1. 
We assumed that H^ ~ 7/8 and 7 ~ 10 using Planckian units. 



There is also an interesting fixed point solution where we can find FiF^ = constant. In this 
case we find that a magnetic field emerges during the inflationary phase due to a logarithmic 
term in the mode expansion of A. Thus the field is driven to be initially large, but decreases 
rapidly as the universe expands. The gauge field then vanishes identically at the end of inflation 
- therefore is unable to act as a seed-field to generate the observed magnetic field today. 
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